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ABSTRACT 

Aims. The super-earth planet GJ 1214b has recently been the focus of several studies, which use the transit spectroscopy technique to 
determine the nature of its atmosphere. Here, we focus on the Ha line as a tool to further restrict the nature of GJ 1214b's atmosphere. 
Methods. We used the Gran Telescopio Canarias (GTC) OSIRIS instrument to acquire narrow-band photometry with tunable filters. 
We were able to observe the primary transit of the super-Earth GJ 1214 in three bandpasses: two centered in the continuum around 
Ha (653.5 nm and 662.0 nm) and one centered at the line core (656.3 nm). We measure the depth of the planetary transit at each 
wavelength interval. 

Results. By fitting analytic models to the measured light curves, we were able to compute the depth of the transit at the three 
bandpasses. Taking the difference in the computed planet-to-star radius ratio between the line and the comparison continuum filters, 
we find A{R p /R ir ) Ha -6S3.5 = (6.60 ± 3.54) x 10~ 3 and &(R p /R*) Ha -662.o = (3.30 ± 3.61) x 10~ 3 . Although the planet radius is found to 
be larger in the Ho- line than in the surrounding continuum, the quality of our observations and the sigma level of the differences (1.8 
and 1.0, respectively) do not allow us to claim an Ha excess in GJ 1214's atmosphere. Further observations will be needed to resolve 
this issue. 
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1. Introduction 

Exoplanetary research is one of the fastest-growing areas of as- 
tronomy. Particularly exciting was the detection of transiting 
planets with a mass lower than Neptune (e.g., Leger et al. 2009 



|Howard et al.|[20TT| i, which has opened up the study of plan 
ets with masses in the range 1.5 to 10 M®, the so-called super- 
Earths. Even more recently, Kepler results have pushed the de- 
tection limit down to Mars-size p lanets (KOI 961.01, 961.02, 
and 961.03; Muirhead et al. 2012). Super-Earths are important 
because they have a mass range that is not present in the solar 
system, they can bring new clues about planetary formation, and 
they are good candidates for future searches for life. 

One of the most-studied super-Earths is GJ 1214b, which 
was discovered by Charbonneau et al. ( 2009 ) as part of a ground- 
based exoplanet search in M dwarf stars (MEarth, |Nutzman &| 



Charbonneau 2008 ). Since this system has a larger planet-to-star 



radius ratio (Rp/R*) than most of the super-Earths detected so 
far (e.g., Corot 7b, Kepler 22b), GJ 1214b has become a very 
interesting candidate for the study of its atmospheric composi 
tion. According to some simulations (Rogers & Seager 2010 



there is a degeneracy in the models that account for the possi- 
ble atmospheric composition of GJ 1214b based on its density 
and irradiation. Thus, the planet can either be composed of a) 
a rocky/ice core surrounded by a primordial atmosphere domi- 
nated by hydrogen and with a relatively large-scale height, b) a 
water and ice core (Waterworld) with an atmosphere dominated 
by water vapor with a small-scale height, or c) GJ 1214b can be 



a rocky planet with an atmosphere formed by outgassing with a 
large atmospheric scale height. 

Transmission spectroscopy has proven to be a successful 
technique to constrain the planetary composition and even detect 
molecules in extra-solar planets (e.g ., |Charbonneau et al.|2 002 ; 
Tinetti et al.|2007 Sing et al.|2011 1. The methodology is based 



on the fact that each element and molecule present in the plan- 
etary atmosphere possesses a different opacity, thus producing a 
slight change in the observed radius of the planet R„ (therefore 
changing the ratio R p /R+) at different wavelengths. By compar- 
ing the measured changes in Rp/R* with the predictions made 
by different atmospheric models, one can infer the composition 
of the planet. 

In the case of GJ 1214b, several studies have been made to 
try to constrain its atmospheric composition through transmis- 
sion spectroscopy (e.g., Bean et al. |2010 Desert et al. 2011a 



|Crossfield et al.|201 1| |Croll et al.|201 l| i,"and disentangle the de 
generacy between the models that account for its density. Using 
spectra taken with VLT/FORS, |Bean et al.| ( |20T0| did not find ev- 
idence of any significant features between 780 nm and 1000 nm. 
This result has been confirmed by Bean et aL](|201 1[ ) and other 
studies (e.g., |Berta et al.|2012[|de Mooij et al.|2012| > that did not 
find evidence of features present in the atmosphere of GJ 1214b. 
This lack of features is interpreted as evidence for the presence 
of a metal-rich atmosphere and/or thick clouds that produce a 
constant absorption across a broad range in wavelength. 

Here, we present a new set of narrow-band Ha photome- 
try of a transit of GJ 1214b on August 17 2011, taken with 
the Optical System for Imaging and low Resolution Integrated 
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Spectroscopy (OSIRIS, Cep a et al.|[2000) i at the 10.4 m Gran 
Telescopio Canarias (GTC), located at Observatorio Roque de 
los Muchachos in La Palma, Canary Islands, Spain. In this study 
of the atmospheric signature of a super-Earth, we use a differ- 
ent approach than in the literature; instead of using broadband 
photometry or low-resolution spectroscopy, we use narrow-band 
imaging with tunable filters. 



2. Data 

2.1. Observations 

The OSIRIS instrument consists of a mosaic of two Marconi 
CCD detectors, each with 2048 x 4096 pixels and a total unvi- 
gnetted field of view of 7.8 x 7.8 arcmin, giving a plate scale 
of 0.127 arcsec/pix. For our observations, we choose the 2x2 
binning mode with a readout speed of 200 kHz (that has a gain 
of 0.95 e-/ADU and a readout noise of 4.5 e-), defining a single 
window of 1.45 x 6.87 arcmin to increase the sampling in the 
sequence. This configuration produced a readout overhead time 
of only 2.9 s. 

We used the OSIRIS Red Tunable Filter (RTF) for our ob- 
servations, which consists of a Fabry-Perot etalon that allows 
for narrow-band imaging, selecting both the central wavelength 
and the width of the filter. In the case of OSIRIS, the user can 
select a filter centered in the wavelength range of 651-934.5 nm, 
with a width of 1.2-2.0 nm. The use of tunable filters (TFs) in 
exoplanetary science has been explored in recent years (e.g., 
Colon et al.pOlO] |Sing et al.|[20TT] |Col6n & Fordl|20TT) . The 
collecting capability of the 10.4 m GTC telescope in combina- 
tion with the flexibility of the TFs makes it possible to obtain 
transit depths at a fixed wavelength with enough accuracy to de- 
tect possible atomic signatures in planet atmospheres. This was 
the case with the first detections of potassium in the extrasolar 
planet XO 2b (Sin g et al.|2011|l and th e study of the same ele- 
ment in HD 80606b ( jColon et al.|2012) . 

As in any Fabry-Perot etalon, each pixel in the CCD in 
OSIRIS RTF images has a slightly different wavelength, increas- 
ing radially from the optical center. This increase thus produces 
"rings" of constant wavelength with respect to the RTF center 
located nearly in the middle of the detector system. As a result, 
stars with different positions will have different observed wave- 
lengths, unless they are at the same distance from the optical 
center. For this reason, the field of view of OSIRIS was rotated 
149.63 deg (North is down-right in Figure [TJ in order to have 
GJ 1214 and one bright reference star located at the same dis- 
tance, about 3.05 arcmin, from the optical center. In this way, 
they can be at the same wavelength, while the rest of the stars 
are at different distances from the optical center and hence at 
different wavelengths. At the distance of 3.05 arcmin from the 
optical center the tunable filter was set up at 658.0 nm, 661.2 nm, 
and 666.9 nm to get the desired wavelengths of 653.5 nm, 656.3 
nm, and 662.0 nm at both the target and reference star, respec- 
tively. For each narrow filter, we used a width of 1 .2 nm and we 
set an exposure time of nine seconds for each band. This gave 
a duty cycle of more than 75%, taking into account the small 
readout time used along the series. The continuum-line tunings 
were cycled during the observation using increasing wavelength 
(i.e., 658.0 nm - 661.2 nm - 666.9 nm). Figure [T]presents one of 
the raw OSIRIS windowed images of GJ 1214 obtained in the 
series, showing the position of the target and reference star used 
in the photometry. 

Data for GJ 1214 were taken in service mode on August 17, 
201 1. Transit observations began at 22:50 UT and ended at 1:00 




Fig. 1. Raw Ha GTC/OSIRIS image. The position of GJ 1214 
and the reference star used in the photometry is shown. 



UT, during which the airmass ranged from 1.27 to 2.37. The ob- 
serving conditions were good, with a clear sky and bright moon. 
Seeing was good and stable (0.7-0.8 arcsec) along the whole 
transit, so that a slight defocus was needed to avoid saturation in 
the reference star. This defocus was produced by changing the 
nominal position for GTC secondary mirror by ±0.1 mm and 
maintained without changes for the three different wavelengths. 
With this procedure, a maximum count level of around 35,000 
- 40,000 ADUs was ensured along the observations at the refer- 
ence star, the brightest star in the field, well below the saturation 
level (65,000 ADUs) for OSIRIS CCDs. At 23:45 UT, near mid 
transit, some guiding problems (caused by telescope defocus) 
resulted in the loss of some images and a slight change in the 
position of the stars (about eight pixels) for the rest of the ob- 
servations. This change, due to the decrease of the radial wave- 
length outwards from the optical center, translates to a difference 
of less than 0.05 nm with respect to the initial wavelength setup 
and thus has no effect in the photometric measurements. The 
RTF also suffers from a calibration effect that depends on the ro- 
tator angle (see OSIRIS user manual However, for the range 
of rotator angles covered along the sequence, the differences in 
the TF tuning for this effect were again lower than 0.05 nm, so 
that no recalibrations were needed during the transit. 

2.2. Data reduction and Photometry 

Photometric reduction was completed using standard proce- 
dures. Bias and flat-field images for each filter were processed 
using the Image Reduction and Analysis Facility (IRAlQ rou- 
tines, while the bias, flat-fielding correction, centroid compu- 
tation, and aperture photometry were performed using specific 
Interactive Data Language (IDL) scripts. 

Finally, only one star was used as reference for the differen- 
tial photometry, since we could only get one bright star at exactly 
the same distance as GJ 1214 from the optical center (that is, at 



1 www.gtc.iac.es/en/pages/instrumentation/osiris.php 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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the same wavelength as GJ 1214). Although tests were made to 
compute the flux ratio bet ween GJ 1214 and other reference stars 
that were available within the OSIRIS field of view (FOV), we 
concluded that our initial reference star was the best option be- 
cause it delivered better results in rms than an ensemble of stars. 

Due to the wavelength variation across the OSIRIS RTF 
FOV, ring-like structures resulting from sky emission lines ap- 
pear in the images. To remove any noise that these structures 
might introduce to the photometry, we used the IRAF package 
which was created to reduce images taken with tunable filter s 
at the Anglo-Australian Telescope TFRecj^] ( |Jones et al. 



2002) 



This artificially dithers each frame in order to create a sky im- 
age, which is then substracted from the original data. The use 
of TFRed, however, did not improve the quality of the photom- 
etry, as noted in previous works (Colon et al. ( 2010[ l), so we 
decided to follow the traditional approach of performing aper- 
ture photometry by using a fixed aperture and a sky substraction 
based on the average counts in an annulus away from the star. 
The aperture and sky ring combination that delivered the lowest 
rms was an aperture of 2.25 arcsec (approximately three times 
the measured FWHM of the objects) and a sky level computed 
using a ring of 2.5 arcsec of width that started five arcsec away 
from the center of the stars. The errors in the photometry were 
computed by the IDL routine APEB^] which uses photon count 
statistics to estimate the error level. However, the computed error 
was smaller than the point-to-point spread observed in the light 
curve, so we adopted the rms for the out-transit data as final error 
estimate. 

After computing the flux ratio between the target and the ref- 
erence star, a second-degree polynomial was fitted to the out-of- 
transit data points in order to correct a trend in the light curve. 
This trend is probably produced by atmospheric extinction af- 
fecting the flux of the stars in a different way due to their distinct 
color. Using the correlation between the flux ratio and the air- 
mass, another polynomial fit was made as a test. However, since 
it delivered similar results in terms of rms of out-transit data and 
curve shape as using the temporal trend did, we decided to use 
the initial correction for the light curve presented here. 

We obtained the universal time of data acquisition for each 
frame by using the recorded headers of the images, which in- 
dicate the opening and closing time of the shutter, in order to 
comp ute the time of mid ex posure. Then using the code writ- 
ten by|Eastman et al. (2010|P] we computed the barycentric ju- 



lian date in barycentric dynamical time (BJD) using the mid- 
exposure time for each of the science images to produce the final 
light curves at the three different wavelengths analyzed here. 



3. Data Analysis 

3.1. Light Curve Fitting 

The light curves obtained for each of the three filters, together 
with the best-fitting synthetic models, are shown in Figure [2] 
As pointed out, the telescope temporarily lost the guiding near 
mid-transit, producing a change in the position of the stars and a 
few outliers in the photometry, which are more noticeable in the 
656.3 nm and 662.0 nm filters. 

To fit the observed light curve to synthetic models, we fixed 
two of the orbital elements, the ratio a/R+ = 14.9749 and the 
inclination of the system i = 88.94° (both taken from Bean 



http ://www.aao . gov, au/local/w w w/j bh/ttf/ a dv_ reduc .html 



4 http://idlastro.gsfc.nasa.gov/ 



http://astroutils.astronomy.ohio-state.edu/time/ 



|et al. |201 1 1 and also assumed an eccentricity of e = 0. Since 
GJ 1214b did not show any evidence of central transit time vari- 
ations ( |Bean et"aL|2010 Berta et al. 2012 1, we fixed the pe- 
riod to P = 1.58040481 days ( |Bean et al.pOll) . The central 
time of the transit (To) was left as a free parameter for each of 
the three filters and so that we could see if there was consider- 
able departure from the predicted To. The predicted central time 
of T Eph = 2455791.49643382 BJ D was calculated u sing a lin- 
ear ephemeris (with Tq taken from Bean et al.|201 \\ . Since our 
data consists of observations taken with very narrow filters, we 
considered that fixing the limb-darkening (LD) coefficient value 
extrapolated from the results found for broadband photometry 
could introduce a source of error in the estimation of the transit 
depth. Thus, we decided to also leave the LD coefficient as a free 
parameter in the fitting process. 

The models of Gimenez (2006) were fitted separately for 
each filter, leaving as free parameters the central time of the tran- 
sit, LD coefficient, and the depth of the transit. Here we have 
chosen to use a linear LD law, instead of a quadratic or nonlin- 
ear coefficients, because the data is noisier than a priori expected 
based on telescope aperture. Thus, we consider that the use of a 
linear LD law is enough to account for the characteristics of the 
curve. 



3.2. Parameter Estimation 

In this section, we describe the bayesian method used to estimate 
the fitted parameters to the light curve. The procedure is similar 
to the one used, and amply detailed, in Berta et al. (2012). 

The first step to estimate the distribution of the free param- 
eters chosen to fit the observed light curve (depth and LD) is to 
minimize the x 2 function that compares the data with the syn- 
thetic models. It is: 



X 



z 



(di - m,f 



cr 



(1) 



where d, is the observed point of the light curve, m,- the point 
taken from the comparison model, and cr the error in the obser- 
vation. Once x 2 1S minimized, we have an initial guess for the 
transit depth and LD coefficient to feed to our bayesian analysis. 

In the bayesian analysis, the Metropolis-Hasting method is 
applied to find the distribution of probabilities of the transit 
depth and LD parameters. The algorithm compares the obser- 
vations (D) with different synthetic models (M), which were 
computed using the transit depth and LD coefficient generated 
randomly with a probability distribution. In our case, and for the 
parameters mentioned before, we adopted a normal random dis- 
tribution centered in the value found by the x 2 minimization as 
well as a standard deviation found by trial and error to see which 
combination of numbers decreased the time of convergence of 
the algorithm. In each step, a comparison is made using Bayes' 
theorem to see if the model in the t + 1-th iteration (M t +i) has 
higher probability than the previous step (M,). This is done by 
evaluating their ratio of probabilities r given by 



P(D\M l+l ) P(M t+l ) 
P(D\M t ) P(M t ) 



(2) 



where P(D\M) is the posterior likelihood of the model, which 
evaluates how good our choice of parameters is compared to the 



data. In our case, we computed the likelihood as in |Berta et al. 
( |20T2l : 



\nP(D\M) = -N\n(scr) 



X~_ 
2s 2 



■ constant , 



(3) 
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where N is the number of data points, x 1 is described by 
Equation [T] and s is a parameter introduced to take into account 
the true noise level in our data. In each iteration, we compute Eq. 
fusing a normal random distribution to generate a new value for 
Tq, s, the transit depth, and the LD coefficient. 

Called the prior probability distribution, P(M) describes the 
knowledge that we have of the space in which our parameters are 
distributed. For the central time (To), depth (k), and s parameter, 
we use Jeffreys prior and for the LD coefficient (c) uniform prior: 



P(M) 



1 



1 



1 



To 

l 



(4) 



where x,„,„, x max with x = [To, k, c, s] represents the limits of 
the domain of the central time, depth, LD, and s respectively. For 
the central time, we use [T Eph -0.00l042,T Eph +0.001042], for k 
[0.008,0.5], for LD [0.1,1.0], and for s [0.5,2.0]. 

Once we have computed the ratio r, we compare it to a num- 
ber taken from a random uniform distribution U in the inter- 
val [0, 1]. If r > U(Q, 1), we accept the present model and set 
M t — M t +\. To ensure that the method converges to the distri- 
bution of probabilities for each fitted parameter, this process is 
repeated several times. In our case, we performed 2.5 x 10 5 iter- 
ations for each light curve. 

Since the algorithm takes time to start sampling the true dis- 
tribution of the parameters, we have to eliminate some of the first 
iterations in a process known as a burn-in period. For the depth, 
central time, LD coefficient, and s, we did not consider the first 
2 x 10 4 iterations. By computing the autocorrelation for each fit- 
ted parameter, we were able to obtain the correlation length of 
the chain and keep only the parameter values that were indepen- 
dent in order to establish their probability distribution. At the 
end, for each parameter we have of the order of 10 4 points to 
compute their distribution. 

Once we obtained the probability distribution of the param- 
eters, we computed the percentiles of the distribution in order 
to obtain the median and standard deviation for each parame- 
ter. The median is the value below which 50% of the observed 
points are found, while the standard deviation was obtained us- 
ing the value by which 68.27% of the data is found (around 
the median). By using this procedure, we are able to find the 
median and standard deviation for a nonsymmetric function, al- 
though most of the fitted parameter distibutions present gaus- 
sian shapes. In Table [T] we give the results of the light curve 
fitting and parameter estimations for the three different filters. 
The values and quoted errors are the median and standard devi- 
ation derived from the probability distribution, respectively. We 
used the median values for the central time, transit depth, and 
LD coefficient to produce the model curves and the s parameter 
to estimate the true level of error by e true = scr (see Figure [2). 
Looking at Table [T] we note that the LD coefficients computed 
have consistent values for all the filters within the relative large 
errors, which they seem to increase to shorter wavelengths. 

The rms of the residuals from the best fit for each light curve 
(Fig. [2) are 3900 ppm (653.5 nm), 3100 ppm (656.3 nm) and 
3200 ppm (662.0 nm). 

As an additional test, we performed a prayer bead analysis 
( |Moutou et al.|2004| |Gillon et aLpOOT] |Desert et al.|201 lb| > in 
each filter (separately) to get another estimation of the fitted pa- 
rameter median values and errors and to see the possible impact 
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Fig. 2. GJ 1214b transit light curves obtained with GTC/OSIRIS 
tunable filters. In the bottom panels, the residuals of the curve fit 
are shown for each filter. 



Table 1. GJ 1214b central transit times, depths, and linear limb- 
darkening coefficients from GTC/OSIRIS data. 



Filter 


R P /R* 


c 




AT (days)" 


653.5 nm 


0.1151 ±0.0025 


0.670 ± 


0.143 


0.00040 ±0.00031 


656.3 nm 


0.1217 ± 0.0025 


0.793 ± 


0.122 


-0.00036 ± 0.00035 


662.0 nm 


0.1184 ±0.0026 


0.917 ± 


0.099 


0.00007 ± 0.00044 



<n) Offset A7o = T Ep it - Tq computed from ephemeris central time 



' Eph 



■■ 2455791.49643382 BJD. 



of systematics in the fitting process. The test consisted of shift- 
ing the residuals from the initial fit (in this case using the values 
from Table [TJ, adding these residuals to the initial fit, and ad- 
justing the resultant light curve again; this process is repeated N 
times, where N is the number of observed data points. 

For all the adjusted parameters, the median values found 
by the prayer bead test were consistent with the ones found by 
the bayesian approach, meaning that these parameters were not 
strongly affected by systematic errors. The errors computed us- 
ing this method were smaller than the ones previously found, so 
we adopted as final error values the ones found by the bayesian 
approach. 
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4. Discussion and Conclusions 



It is easily seen in Figure [2] that the light curves show a fair de- 
gree of variability throughout the observations. Given the good 
quality of the weather and sky conditions during the observa- 
tions, this noise may have arisen from an imperfect correction 
of atmospheric effects due to the use of only one reference star 
located relatively far from the target (~ 5.5 arcmin) and the high 
airmass (> 2) at the end of the observations. This photomet- 
ric dispersion could also have been produced by stellar activity; 
GJ 1214 has previously shown evidence of some level of activ- 
ity in the form of spots and low-energy flares ( |Kundurthy et al. 



2011 



et al. 



, although it was first classified as an inactive star (Hawley 
1996 ). This phenomenon could have affected our photom- 



etry since we focused on the Ha line, which has been used as an 
indicator of stellar activity. The same effect could be present in 
the reference star used in the differential photometry. 

In Fig. 3] we compare our results with previous works on 
GJ 1214b. It shows how well this planet has been studied since 
its discovery, and it possesses one of the most complete trans- 
mission spectra to date. Our data presents a bigger scatter when 
compared with the work of Bean et al. (201 1 1, which gave op- 



tical and near-infrared data taken at several facilities. The opti- 
cal observations were made with the VLT FORS spectrograph, 
where the spectrum was binned into 20 nm wide wavelength in- 
tervals (a factor of 16.66 coarser spectral resolution than our RTF 
data), spanning 610 nm to 850 nm. The interval that contains our 
three data points (650 nm - 670 nm) presents a planet-to-star ra- 
tio of R p /R+ = 0.1167 + 0.0011, which differs from our value 
for the Ha line of Rp/R* = 0.1217 + 0.0025. However, if we 
consider our two RTF setups in the continuum near the Ha line, 
the retrieved planet-to-star radius ratio is fully consistent with 
the value found by Bean et al. ( 2011| l within the error bars. 

If we take the difference in the computed planet-to-star ra- 
dius ratio between the line and the comparison continuum fil- 
ters, we find that A(^/^) Hff -653.5 = (6.60 ± 3.54) x 10~ 3 and 
A(tf / ,//?*)/ /u ,_662.o = (3.30 ± 3.61) x 10~ 3 . The comparison be- 
tween Ha while the blue continuum (653.5 nm) is consistent 
at 1.8cr and between Ha and the red continuum (662.0 nm) is 
consistent at 0.9<x level. As a result, we did not find a statisti- 
cally significant detection of Ha. Still, the larger-than-expected 
planet radius in the Ha line probably deserves new observations, 
aiming at higher SNR, to finally solve this issue. 

This nondetection of Ha is coherent with previous lower res- 
olution results that do not find evidence of features presented in 
the atmosphere of GJ 1214b. In addition, our results can be used 
to constrain limits on the transmission spectrum GJ 1214b at 
high resolution. 
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Fig. 3. GJ 1214b transmission spectra obtained with 
GTC/OSIRIS tunable filters compared with other works. 
Top panel: full range of observed wavelengths by different 
authors, Bottom panel: optical observations (600 nm - 700 nm); 
the horizontal bars represent the width of the filters used in each 
observation. 
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